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SUMMARY 

This paper considers the reported attack ratio arising from outbreaks of influenza in enclosed 
societies. These societies are isolated from the wider community and have greater opportunities 
for contact between members which would aid the spread of disease. While the particular kind of 
society (prison, care home, school, barracks, etc.) was not a significant factor in an adjusted 
model of attack ratio, a person's occupation within the society was. In particular, children and 
military personnel suffer a greater attack ratio than other occupational types (staff, prisoners, 
etc.). There was no temporal trend in final attack ratio nor, with the exception of 1918, do 
pandemic years show abnormal attack ratios. We also observed that as community size increases, 
the attack ratio undergoes steep nonlinear decline. This statistical analysis draws attention to how 
the organization of such societies, their size and the occupations of individuals within them affect 
the final attack ratio. 

Key words: Analysis of data, community epidemics, infectious disease control, influenza, mathematical 
modelling. 



INTRODUCTION 

The transmission of respiratory infectious diseases 
and, a usual epidemiological corollary of this, the 
final clinical attack ratio at the end of epidemics is 
known to be influenced by factors such as population 
density and the patterns and intensity of population 
mixing at different spatial and social scales [1-3]. In 
the case of influenza, secondary attack ratios in house- 
holds have often been shown to decrease with increas- 
ing household size [4] and attack ratios are also 
generally expected to be higher in schools than in 
the wider community [5], Here we use the dimen- 
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sionless terminology attack 'ratio' to represent the pro- 
portion of the population who fall ill with influenza 
over the whole reported outbreak. 

Enclosed societies form an important spatial scale 
between the household and the general population. 
These societies are diverse both in terms of their func- 
tion and the extent to which social interactions are 
encouraged or allowed between members and sub- 
groups. Some form part of an institution (prisons, 
schools, the military, etc); others are isolated by natu- 
ral factors. Individuals within these societies may be 
sufficiently close that repeated contact with the same 
people each day may enhance the mechanisms for dis- 
ease transmission. In certain societies the importation 
of disease may have a significant impact on the mem- 
bers of the society. Following importation, diseases 
that are subclinical or difficult to differentiate on 
clinical grounds may spread quickly before local 

t subject to the conditions of the Creative Commons 
ises/by-nc-sa/3.0/>. The written permission of 



108 T. J. R. Finnie and others 



Table 1. Descriptive statistics of available risk factors for an enclosed 
society attack ratio 





Number of 


Number of 








institutional 


individuals 


Attack 


Median 


Characteristic 


groups (n - 83) 


(w = 39602) 


ratio (%) 


(IQR) 



Year 
1890 
1918 
1920 
1957 
1972 
1976 
1978 
1980 
1983 
1984 
1986 



1997 
2000 
2001 
2002 
2004 
2009 

Institution category 
School 
Barracks 
Children's home 
Hospital and care home 
Isolated community 
Prison 
Ship 

Occupation type 
Children 

Military personnel 
Patients and residents 
Prisoners 
Staff 

Group size 



37 
10 
1 
3 
2 
1 
2 
1 
1 
1 
1 
2 
5 
2 
2 
1 
2 
9 

24 
3 
9 

19 
1 

18 

9 

30 
9 
12 
13 
19 



8121 
9353 
1000 
223 
4148 
859 
1107 
859 
859 
555 
859 
1804 
211 
41 
80 
3000 
1481 
5042 

9115 

8239 
1366 
4330 
3200 
9428 
3924 

14120 
11560 
3279 
9084 
1559 



20 
32 
44 
73 
22 
24 
52 
23 
29 
14 
21 
44 
40 
22 
31 
3 
17 
17 

36 
18 
27 
16 
18 
13 
53 

32 
28 
15 
13 
19 



185 (49-610) 



IQR, Interquartile range. 



behaviour is changed to mitigate the outbreak [6-10]. 
Furthermore, some of these societies may contain 
members that are particularly susceptible to infection 
or are severely compromised by infection (immuno- 
compromised people). 

The aim of this study is twofold, first, to analyse the 
reported outbreaks of influenza in enclosed societies 
over the past 120 years to determine the quantitative 
evidence for the hypotheses that attack ratios can be 
higher in certain types of enclosed societies than 
others; and second, to test the further hypothesis 
that in enclosed societies with larger populations 
there is a tendency for lower final attack ratios. It 



was recognized that there would be confounding fac- 
tors to such an analysis, as, for example, the strain 
of influenza implicated in the outbreak and the past 
history of influenza in a community. These factors 
are unreported in most published data. 

METHODS 
Data 

The data collected and described by Finnie et al. [11] 
provide information on influenza outbreaks in 43 
enclosed institutions over the years 1890-2009 
(Table 1). The outbreaks are distinguished by the 
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Table 2. Estimates of fixed effects for an enclosed society attack ratio. All regression models have year and 
institution as random effects. Referent is the subgroup of a categorical variable against which other subgroups 
are compared in the regression mode 



Characteristic 


Unadjusted OR 
(95% CI) 


P value 


Adjusted OR 
(95% CI) 


P value 


Institution category 
School 
Barracks 
Children's home 

nuapilal ailU Calc IlOIIlc 

Isolated community 

Prison 

Ship 


Referent 

0-17 (0-05-0-59) 

0-55 (0-17-1-77) 

0-38 (0-06-2-55) 
0-24 (0-11-0-55) 
0-89 (0-39-2-01) 


0-02 


Referent 

0-67 (0-12-3-73) 

0- 58 (0-20-1-65) 

9-18 ffl-R^ S-7T> 

Z lu 03 — J 1 Jf 

1- 49 (0-27-8-24) 
1-06 (0-44-2-57) 
1-44 (0-44-4-70) 


0-40 


Occupation type 
Children 

Military personnel 
Patients and residents 
Prisoners 
Staff 


Referent 
0-56 (0-25-1-29) 
0-36 (0-23-0-54) 
0-35 (0-22-0-54) 
0-47 (0-33-0-68) 


<0-01 


Referent 
0-73 (0-19-2-86) 
0-12 (0-07-0-20) 
0-26 (0-16-0-44) 
0-16 (0-10-0-25) 


<0-01 


Group size (In) 


0-80 (0-75-0-86) 


<0-01 


0-63 (0-57-0-69) 


<001 



OR, Odds ratio; CI, confidence interval. 



occupational group of the institution's population 
affected and the size of this population, as well as 
the final number of individuals who acquired in- 
fluenza. The data summarize 9737 cases of influenza 
occurring in a total population of about 39600. 

Institutions with regimens that permit more fre- 
quent, longer term and more 'intimate' contacts 
between individuals would generally be expected to 
promote increased overall attack ratios. The contact 
regimen may be related to the purpose of an insti- 
tution and institutions were accordingly classified 
into one of seven categories on the basis of their func- 
tion (Table 1). In addition, as, a priori, we would 
expect there to be differences in attack ratios due to 
different activities, so each of the population groups 
within an institution was categorized into a broad 
occupational type (Table 1). 

Statistical analysis 

Attack ratios are correlated within years and insti- 
tutions and we accordingly fitted mixed-effects logistic 
regression models [12]. Year and institution were 
modelled as normally distributed random effects 
with a mean of zero. Year was effectively a proxy 
for influenza strain. Pandemic years might be expected 
to have larger outbreaks due to larger susceptible 
groups in wider populations and so the normality 



assumption for year was tested using the Anderson- 
Darling test (a standard test for a sample of data 
being derived from a normally distributed population) 
[13]. Institution was a proxy for local response and 
other similar factors. Group size, occupational type 
and institution category were included as fixed effects, 
at first in separate models and then in a joint multiple 
regression to obtain adjusted estimates of their effects. 
Occupation type was to some extent confounded with 
institution category and we were not able to examine 
the interaction between these two effects due to data 
limitations. Group size follows a right-skewed distri- 
bution and was allowed to vary on a natural logarith- 
mic scale in the model as initial investigation revealed 
this provided a better fit to the data than a linear 
term. Attack ratio data were converted to a binomial 
response variable and analysed with the lme4 package 
[14] in R using a logit link function. The significance 
of the random and fixed effects was assessed using 
likelihood ratio tests. 

Table 1 shows that data are included for one iso- 
lated community consisting of 3200 individuals. 
Data on occupational groups within this population 
were not available and so for the purposes of statisti- 
cal analysis the entire community population was 
assigned to the 'children' occupation category. The 
sensitivity of results to this assignment were investi- 
gated by omitting the community from this category 
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Fig. 1. Attack ratio and group size by broad occupational type. Estimates of effect sizes drawn from multiple regression 
models using mean value of zero for random effects, two occupational categories as shown, and reference institution 
category (school). 



and placing it in an occupational category of its 
own, for the single fixed-effect regression, or omitting 
it from the analysis altogether for the multiple 
regression. Resultant changes in model coefficients 
and significance levels were then examined. 

RESULTS 

Unadjusted estimates from the single covariate mixed 
effects regressions of institution category, occupation 
category and group size on attack ratio are given in 
the second column of Table 2. Adjusted estimates of 
these effects from the multiple regression are given 
in the fourth column of Table 2. Although all fixed 
effects are significant at the 5% level when considered 
singly, only occupation category and group size 
remain significant at the 5% level in the multiple 
regression. The random effects of year and institution 
were significant at the 5% level in all models. The 
normal distribution was found to be a reasonable 
approximation to the year random effect as the 
Anderson-Darling test did not reject the hypothesis 
of normality at the 5% level. 

The unadjusted odds ratios show that the odds 
of suffering from influenza in an enclosed institution 
are on average lower for individuals in barracks, 



children's homes, hospitals, isolated communities, 
prisons and ships than they are for individuals in 
schools (P = 0-02) (Table 2). The estimates for 
barracks (OR 0-17, 95% CI 0-05-0-59), hospitals and 
care homes (OR 0-34, 95% CI 0-14-0-86) and prisons 
(OR 0-24, 95% CI 0-1 1-0-55) are significantly different 
from the baseline estimate for schools. Occupation 
type within an enclosed institution is, however, more 
significant than institution category per se as an unad- 
justed fixed effect (_P<0-01). The odds of acquiring 
influenza are on average much lower for patients 
and residents (OR 0-36, 95% CI 0-23-0-54), prisoners 
(OR 0-35, 95% CI 0-22-0-68) and staff (OR 0-47, 95% 
CI 0-33-0-68) than they are for children. These esti- 
mates are not sensitive to the occupation type assigned 
to the isolated community as the estimates of effect 
size vary by an average of only 4% between the 
alternative assignments (log scale). The unadjusted 
effect of group size also indicates that the odds of 
acquiring influenza decrease by 20% for every one 
unit increase in the log of group size (OR 0-80, 95% 
CI 0-75-0-86). 

Occupational type and group size remain significant 
in the multiple regression model of attack ratios 
(i , <0 01) but institutional category is not significant 
(,P = 0-40). Patients and residents in enclosed 
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Fig. 2. Magnitude, direction and confidence of effect of 
year, modelled as a random variable, on the attack ratio. 
Estimates taken from an adjusted model. Points show best 
estimate of the effect with lines showing the 95% con- 
fidence intervals. Years where the effect is positive had a 
greater attack ratio than would be expected given the 
other modelled effects, while those that are negative had a 
smaller attack ratio than would be expected. 



institutions have odds of acquiring influenza on aver- 
age only 12% of those of a child in an enclosed insti- 
tution (adjusted OR 012, 95% CI 0-07-0-20). Similar 
odds ratios are seen for prisoners (adjusted OR 0-26, 
95% CI 0-16-0-44) and staff (OR 0-16, 95% CI 0-1- 
0-25). These estimates are not sensitive to the inclusion 
of the isolated community observation as the coeffi- 
cient estimates varied by an average of only 0-1% 
between the observation omitted model and the 
model in which isolated community was allocated to 
the 'children' occupational type. The odds ratio for 
group size in the adjusted model falls to 0-63 (95% 
CI 0-57-0-69), indicating that after adjustment for 
occupation and institution category the odds of 
acquiring influenza in an enclosed institution decrease 
by 37% for every one unit increase in the log of group 
size (OR 0-63, 95% CI 0-57-0-69). The combined 
effects of group size and occupational type are demon- 
strated in Figure 1 which gives the estimated expected 



decrease in attack ratio associated with an increase in 
community size for two broad categories of occu- 
pational type. These two categories were created by 
amalgamating those occupational types with statisti- 
cally similar effect sizes (statistical characteristics). 
Figure 1 shows that for a given group size the occu- 
pational types of children and military personnel are 
on average associated with higher attack ratios than 
the occupational types of patients, prisoners, and 
staff. Figure 1 also illustrates that as community size 
increases, the attack ratio undergoes a steep nonlinear 
decline. 

Estimates of the random effects obtained for year 
and institution are given in Figures 2 and 3, respect- 
ively. Year acts as a proxy for both the strain of, 
and local susceptibility to, the influenza virus and 
the methods of outbreak control deployed. Pandemic 
years might a priori be expected to have higher attack 
ratios due to higher susceptibility but we found no 
clear temporal trend in the random effect observed 
by year (Fig. 2). The institution random effect will 
reflect the general nature of the institution's regimen, 
including local response to outbreaks. Groups within 
more relaxed regimens suffer greater attack ratios 
given similar institution type, occupation type and 
size (Fig. 3). 

DISCUSSION 

Review of past outbreaks illustrates the difficulty in 
determining the attack ratio of influenza in enclosed 
societies. Each outbreak is essentially the result of a 
series of random events marked by the choices and 
contacts made by the people involved, any prior 
immunity, the response of public health officials and 
the susceptibility and size of the population into 
which the disease is brought. In order to fully under- 
stand the enhanced attack ratios in context would 
require contemporary measures of the wider commu- 
nity attack ratio, which is very difficult to source due 
to the nature of reporting behaviour. 

Notwithstanding these difficulties, from historical 
data we have examined the hypotheses that attack 
ratios can be higher in certain types of enclosed 
societies than others, and that larger populations in 
enclosed societies are associated with lower final 
attack ratios. Our results indicate it is occupation 
type within an institution, rather than institutional 
category per se, that is a significant driver of final 
attack ratio. We also found that attack ratio decreases 
as enclosed community size increases. 
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Fig. 3. Magnitude, direction and confidence of effect of individual institution on the attack ratio modelled as a random 
variable. Estimates taken from an adjusted model. Points show best estimate of the effect with lines showing the 95% 
confidence intervals. Institutions where the effect is positive had a greater attack ratio than would be expected given the 
other modelled effects, while those that are negative had a smaller attack ratio than would be expected. 



The significance of occupational type as a predictor 
of attack ratio may be seen as a reflection of the 'strict- 
ness' of the regimen across various institutions. Child- 
ren and military personnel form an occupational type 
that is subject to fewer restrictions and/or exhibits 
more mixing on average than the occupational type 
categories of patients, prisoners and staff. The latter 
three categories have a lower attack ratio, on average, 
than the former, visualized as the difference between 
the lines in Figure 1 . The effects of community size 
and regimen strictness are different facets of the 
same underlying process, the ability of people within 
the community to make contacts and have inter- 
actions which spread the disease. Additionally, larger 
communities may have logistical challenges in diagno- 
sis of members that are suffering illness, particularly if 
they are self-reporting. 



We found no clear trend over time in the random 
effect observed by year (Fig. 2). Attack ratio would 
not necessarily be expected to reflect advances in 
medical technology, unless antiviral counter-measures 
were in widespread use during the outbreaks. Had 
mortality ratio been used then a declining trend 
could have been expected but there were insufficient 
mortality data to perform this analysis. 

The particular nature of an institution's regimen 
also has a large part to play in determining the final 
attack ratio. Ceteris paribus, institutions where there 
is greater interaction have many times the attack 
ratio of similar organizations which are more restric- 
tive in their mixing patterns. For example the largest, 
positive, random effect for institution of 1-33 (Fig. 3) 
gives a point estimate for attack ratio of 79% for an 
observation in the reference categories with median 
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group size of 185 and a value of zero for the year ran- 
dom effect. The same observation with the minimum 
estimated random effect for institution of —1-9 has a 
point estimate for attack ratio of only 13%, a differ- 
ence in attack ratio potentially caused by an insti- 
tutional regimen of 66%. 

The reports analysed here of influenza in enclosed 
societies show two clear trends in attack ratio. The 
first is as the size of a community increases, the attack 
ratio decreases. This is likely to be caused by scale 
effects within the community. In the extreme case of 
a very large community we would expect the attack 
ratio to approach that of the general population 
whereas in these small communities it can be signifi- 
cantly greater. The second trend is that attack ratio 
depends more on the occupation of the individuals 
within a community than on the type of community 
itself, so children and soldiers are at greater risk of 
attack than those in other reported occupations 
even within the same society. This is likely to be an 
effect of the differing nature of the contacts these 
groups make. Overlaid on these broad trends is the 
influence of the regimen of a society. Where broadly 
similar societies are attacked those with regimens 
that are more restrictive in allowing contacts are 
likely to suffer fewer infections than those that are 
more open. 

Perhaps surprisingly there was no discernible tem- 
poral trend in the random effects for year in which 
an outbreak took place nor were there markedly 
higher random effects for pandemic years, with the 
exception of 1918. Pandemic years are spread through 
the range of attack rates. Attacks in 1918 appear 
significantly worse but otherwise, attacks in enclosed 
societies during pandemics are no worse than at 
other times. There are more reports from pandemic 
years and probably more societies are attacked but 
inside an enclosed society the attack ratio for a pan- 
demic strain seems no different to other strains. This 
might be anticipated as an enclosed society may well 
be as naive to a seasonal strain, due to rare importa- 
tions, as society as a whole is to pandemic strains. 
More detailed analysis of this effect requires measures 
of the attack ratios both within the enclosed society 
and the wider community. In the future, serological 
evidence of random cross-sections of the population 
may be compared to serological evidence from local 
outbreaks [15]. The collection and reporting of this 
information would assist in quantifying both relative 
attack ratios and the risk of disease importation to 
these enclosed societies. 
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